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Dynamic correlation and response functions of classical and quantum systems in thermal equi- 
librium are connected by fluctuation-dissipation theorems, which allow an alternative definition of 
their (unique) temperature. Motivated by this fundamental property, we revisit the issue of ther- 
malization of closed many-body quantum systems long after a sudden quench, focussing on the 
non-equilibrium dynamics of the Ising chain in a critical transverse field. We show the emergence 
of distinct observable-dependent effective temperatures, which rule out Gibbs thermalization in a 
strict sense but might still have a thermodynamic meaning. 
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Introduction. The development of experimental tech- 
niques which prevent dissipation in quantum many-body 
systems has triggered increasing interest in the non- 
equilibrium dynamics of such closed systems. The uni- 
tary non-equilibrium dynamics of a system initially pre- 
pared in a state which is not an eigenstate of its Hamil- 
tonian is called a quantum quench. Basic questions as to 
whether a stationary state is reached and how this can 
be characterized naturally arise. These questions have 
been addressed in a number of simple models, includ- 
ing the one-dimensional systems reviewed in Refs. [l|, y]- 
Early studies led to the following picture: Non-integrable 
systems should eventually reach a thermal stationary 
state characterized by a Gibbs distribution with a sin- 
gle temperature. Integrable systems, instead, are not 
expected to thermalize but their asymptotic stationary 
state should nonetheless be described by a so-called gen- 
eralized Gibbs ensemble (GGE) with one effective tem- 
perature for each conserved quantity [3|-l6|. Interestingly 
enough, depending on the specific quantity and the sys- 
tem's parameters a Gibbs ensemble turns out to capture 
anyhow some relevant features of the non-equilibrium 
dynamics of integrable systems [7[. In particular, ob- 
servables that are non-local in the quasi-particles display 
numerically the same relaxation scales as in equilibrium 
with a suitable effective temperature, at least for small 
quenches [g, |7[. Local quantities instead do not, with 
possible exceptions for quenches at criticality. 

Our purpose is to revisit the debated issue of thermal- 
ization in closed quantum systems with tools developed 
for the study of classical and quantum dissipative glassy 
systems. The analysis of thermalization in closed quan- 
tum systems focused so far on the property that expec- 
tation values of quantities — such as (i) the conserved 
energy and two-point correlation functions depending on 
either (ii) one or (iii) two times — should behave, at 
long times, as the corresponding averages calculated on 
suitable statistical ensembles. However, an equally im- 



portant property of thermal equilibrium states of both 
classical and quantum systems is the validity of model- 
independent fluctuation-dissipation theorems (FDTs) [8| , 
which relate linear response and correlation functions in- 
dependently of their functional form. Focusing on two- 
time quantities, we investigate thermalization issues from 
this perspective. 

Before getting into the technical details, let us explain 
why fluctuation-dissipation relations (FDRs) should be 
more relevant to thermalization issues than the precise 
functional decay of observables. Take phase separat- 
ing systems as an example. The expectation values of 
one-time quantities - such as the energy density (i) or 
observables of the type (ii) - reach equilibrium values, 
suggesting the equilibration of the sample at the bath 
temperature. This is, however, incorrect as proven by 
the fact that observables of type (iii), such as delayed 
density correlations and linear responses, decay in more 
than one dynamic scale, and in the slowest they do al- 
gebraically, as opposed to the typical exponential equi- 
librium relaxation. In spite of this, one can still define a 
bona fide effective temperature from FDRs that link cor- 
relations to their associated linear response 19| as long as 
one distinguishes different pairs of observables and the 
time-scales in which they evolve. Indeed, some observ- 
ables basically ignore others — they do not interact — 
equilibrate quickly with the environment and are char- 
acterized by its temperature. This is the case of par- 
ticle velocities. Positions, instead, do not equilibrate 
with the bath but acquire the same effective tempera- 
ture, i.e., they partially equilibrate, in their own com- 
mon regime of relaxation. This notion applies to many 
other systems with slow dynamics and constitutes the 
basis of, e.g., a consistent thermodynamic picture of non- 
equilibrium glassy dynamics as realized in mean-field the- 
ory [lO| . A quantum finite-dimensional example with ef- 
fective temperatures are electronic tight-binding dissipa- 
tivc one-dimensional models driven by external fields \\A 



We explore up to which extent these ideas carry over to 
the non-equihbrium dynamics of a highly excited closed 
quantum system. For concreteness, we focus on criti- 
cal quenches in which the system's parameters for times 
t > are tuned to be at an equilibrium (quantum) criti- 
cal point. The detailed analysis of the effective tempera- 
ture based on FDRs in classical thermal quenches [lOl |l2| 
demonstrated that such a parameter has a thermody- 
namic meaning for very late epochs only [l3|, when one 
of the involved times is much longer than the other (see, 



14 



however, Ref. 

The model. To illustrate our approach, we focus on 
the simplest quantum integrable interacting model, the 
transverse field Ising chain [l| 



Hr 
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H+i 



Taf 



(1) 



with periodic boundary conditions and even L. The 
Pauli matrices a^'^ satisfy the SU(2) algebra on the 
same site i and commute on different sites. In what fol- 
lows we set J, h,kB ~ 1 and wc measure time in units 
of h/J and the temperature T in units of J/ks- The 
model is exactly solved by a Jordan- Wigner transforma- 
tion to free fcrmions followed by a Bogoliubov rotation 
in momentum space |15| . The energy of the elemen- 
tary fcrmionic excitations with momentum k is efc(r) = 
2[r2 - 2Tcosk+ l]i/2. For r = and L ^ oo a quan- 
tum critical point at F = 1 separates a paramagnetic 
phase (PM, r > 1) with {af) = 0, from a ferromagnetic 
phase (FM, F < 1) with spontaneous symmetry break- 
ing (erf) ^ and long-range order along the x direction. 
(af) ^ for all F > 0. 

The system is prepared at i = in the ground state 
\ipo) of Hto: while it subsequently evolves with 'Hr=i, 
i.e., at the critical point. The quench from Fg to F 
injects an extensive amount of energy into the system 
which is henceforth conserved. After a transient (stud- 
ied in Refs. 16l |l7| for the chain with free boundaries) 
the system reaches an asymptotic stationary regime. A 
crucial quantity in the description of the dynamics is the 
difference A^ between the Bogoliubov angles diagonaliz- 
ing Hr and Hto ■ 



... A rr r ^ 4[rro-(F + Fo)cosfc + l] 

cosAfc(F,Fo) = 

efe(i )efc(io) 



(2) 



Afe encodes the dependence on the initial state and fixes 
the non-thermal statistics of the excitations created at 
t = 0. 

A criterion that has been used to define an effective 
temperature T^{r, Fq) is to require that the energy after 
the quench — quantity (i) in the Introduction — equals 
the average over a thermal equilibrium ensemble [l|, |7| . 
For model ([T]) this implies [7[ 
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which results in the T^ shown in Figs. [T] and [3] with a 
dashed black line. Requiring, instead, the integrand in 
Eq. ^ to vanish defines the mode-dependent T^g{r, Fq) 
oftheGGE 1,0. 

FDT. We focus on the symmetrized and anti- 
symmetrized two-time correlations of two operators A 
and B in the Heisenberg representation, Anit) = 

^i'Hrtj^Q—i'Hrt 



C^^{t + to,to) 



{'iPo\[AHit + to),BHit„)]±\^o), (4) 



where [X,Y]± = {XY ±YX)/2. More precisely, we 
consider connected correlations that we still denote by 
C^^ . While C^^ approaches the classical correlation 
function for /i —> 0, C^^ is related to the linear response 
function R^^ through the Kubo formula R^^{ti,ta) = 

5{AH-h(t)B{ti))Hl5h{to)\^^^= {2i/h)C^^{ti,to)e{h - 

to) [S| which is valid in and out of equilibrium [9{t < 0) = 
and 6{t > 0) = 1]. The asymptotic stationary regime 
is formally defined by the limit to ^ oo. When it is 
physically established depends on the observable and Fq. 
Natural choices for A and B arc the order parameter af 
and the transverse magnetization a^ . The af autocorre- 
lation functions, Ci^{t), arc non-local with respect to the 
quasi-particles while those of af, C|.(i), are local in the 
same variables [l|. 

In Gibbs equilibrium at inverse temperature (3 the 
FDT connects C^^ and R'^^ , via the model-independent 
time-domain and frequency-domain relations [Sl 






iaTL\v{pkuj/2)C^^{Lo), (5) 



h IniR-'^^ (w) = tanh(/3fta;/2) C^^ [u) , (6) 

respectively. We reinstated h to make the classical limit, 
E^^(f) = -/3dC^^(i)/di e{t), transparent. The FDR 
definition of an effective temperature amounts to replac- 
ing /3 by /3eff(a;) in Eqs. (O and ^ [with a fo-dependence 
in non-stationary (glassy) cases], i.e., 

h ImR^^iu) = tanhiP^i {Lu)hLu/2] C^^ {uj) . (7) 



This is the definition that we shall repeatedly use below. 
Before presenting our results let us summarize what is 
known about C^'^ for F = 1. In equilibrium (Fq = F) 
(V'olo'f (t + to)<^iito)\^o) decays algebraically as |t|^'^/^ 
at T == and as |t|-i at finite T 0. Out of equilib- 
rium the decay of C^ is |t|~^ [l6|. For the special case 
of a fully polarized initial condition (Fq = oo) CI fol- 
lows the same \t\~'^ decay, as can be inferred from the 
results in Ref. [2[. Instead, a generic exponential re- 
laxation of (i/jolcrf (io)o'j (to)|'0o) was argued in [18[ us- 
ing semi-classical methods and later shown to hold ex- 
actly [6|. This is in contrast to the power- law decay of 
the T = equilibrium order-parameter spatio-temporal 
correlations. As far as we know, Cf_ has not been ana- 
lyzed so far. Here we complete this picture by calculating 



C^ and Cj_ for generic To. Wc also study C^ , where 
AB = M = j:l,a!/L. 

The transverse local magnetization. For Fq 7^ F = 1 we 
found that C^ and R^ decay as a sum of power laws of t 
and i + 2io- Thus, no characteristic time can be identified 
and, in addition, one cannot compare these functions to 
the thermal ones in order to define an effective tempera- 
ture, as done in Refs. |30,ll3 ^^^ other observables [type 
(iii) of the Introduction] . Taking to -^ 00 one finds the 
stationary relaxation 

C^(t) = -(87r<2)-i cos 8t + 0(^-3) ^ (8) 

R'{t) = (4^t2)"i [T-i - sinSf] + 0(t-3) , (9) 

with T = [(1 + Fo)/(l - Fo)]^. The complementary 
analysis in the frequency domain allows us to define a 
frequency-dependent T^f^ via the FDR in Eq. ([7]). The 
function T^g{uj) is shown in Fig. [T] for Fq = 0.3 as a 
(red) solid line and it has to be compared with the con- 
stant value obtained from Eq. ^, as in Ref. [7[, shown 
as a black dashed line. The asymptotic regime corre- 
sponds to the limit w — >■ 0, zoomed in the inset, in 
which /3^g = l/T^g diverges logarithmically with the law 
-l/2(l-l/T)/[H-(T-2)arctan(\/T~^)/VT"^] Inw 
(green dashed line). We conclude that, although (erf) 
takes a thermal value [7[ , the dynamics of af is not com- 
patible with an equilibrium thermal behavior. For other 
values of Fq, T^g still vanishes at w = and uj = w^ax = 8 
but it is not concave for Fq > 0.35 [19|. For increasingly 
narrower quenches with Fq — > 1, (3^^ — > c« uniformly 
over all frequencies. 

The transverse global magnetization. The long-time sta- 
tionary decay of the global magnetization correlations are 
even slower than the one of C|. (t) [20| : 

Cl\t) = (8^At3/2)-l gi^(g^ _ ^/4) ^ 0(i-5/2)^ (^Q) 

R^\t) ^ -(4^/^t3/2)-l cos(8i - 7r/4) + 0{t~^l'^). (11) 

The leading-order decay i~^/^ is the same as in equilib- 
rium at finite T. However, while the prefactor depends 
upon T in equilibrium |2l| , out of equilibrium the depen- 
dence on Fq appears only at the next-to-leading order, 
i.e., the long-t limit of C^(t) does not retain memory of 
the initial condition. 

The (blue) dash-dotted line in Fig. [1] is T^g{uj) as 
obtained from the FDR dZ]), applied to Cf. T^^{u) 
approaches 2l\jt at low frequencies, it vanishes for 
ijj = Wniax and becomes non-monotonic for Fq <^ 0.27, 
developing a shallow local maximum. Naively, one may 
expect to recover this value by treating the time-domain 
FDR in the long-i limit as follows. Replacing /3 by a 
constant effective value /3*g in the rhs of Eq. ([5]) , the in- 
tegral can be written as series of odd time derivatives of 
C^it). Inserting Eqs. dTO]) and ((TTl) in the rhs and Ihs 
of this expression, respectively, yields 1 = tanh(4/3*ff) T^ 



Tjf(a; ^ 0) 7^ Tig indicates that ,3^^|(w) cannot be ap- 
proximated by an average constant in the integral. In- 
deed, since only the derivatives of the oscillating factor 
in PH)) contribute to the leading order of Eq. ([S]), T^g 
is the one of the oscillatory frequency, which coincides 
with the threshold value Wmax and, for w — > cjnmx, /^^^ 
diverges as /3pg(w) — — ln(wniax — w)/4. Such threshold 
results from the maximum of the dispersion relation and 
the quadratic dependence of M on the fermionic excita- 
tions. It is therefore unclear whether T^g({xi = W) can be 
recovered from the FDR in the time domain. 

Interestingly enough, it turns out [19| that each fre- 
quency w selects a mode fc such that 2efe = w and T^g (w) 
defined here from the FDR ([7]) coincides with the tem- 
perature T% of the GGE. 
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FIG. 1: (Color online.) Comparison between the effective 
temperature defined via the frequency-domain FDR applied 
to erf (red solid line) and M (blue dash-dotted line) , and the 
energy definition (black dashed horizontal line) for Fq = 0.3. 
The inset highlights the logarithmic divergence of Tgjj(cj) for 
oj -i>0. 



The order 'parameter. Equations ©, dH), dTO]) and pi|) 
are invariant under Fq i— >■ Fq because cosA/j(l,Fo) = 
cos Afe(l,F,^"^) [see Eq. ([2])] is the sole quantity bringing 
about the dependence on Fq in the stationary limit to -^ 



00 of C: 



M-. 



In the same limit and for F ~ 1 we 



find numerically that this invariance also holds for C^. 
Henceforth we restrict to quenches originating from the 
FM phase. We computed C\ for a chain with L = 10'^ 
and to = 10 with the methods employed in Refs. [71, ll5|, 
Il6l | . Our numerical results are fitted very accurately by 

Cl{t)^e-iAc[l + act--2 sin(4i + 0)], (12) 

i?^(f) ~e--^fl[l-afl i-5 cos{At + (j))], (13) 

with (numerically) the same rate t~^ as the one t~^ = 
— /J^(dfc/7r)[defe(F)/dfc] Incos Afe analytically proved [6| 
to characterize the exponential long-time decay of 
{il,'o\iTf{t)aJ{t)\^/jo). The expression for r finds fur- 
ther support from the fact that with the substitution 
cosAi, — >■ tanh(/3efe/2) one recovers the equilibrium 
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c» and therefore T*fr 



0. The fact that 
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23[ in this as well as other statistical aver- 



Although several fitting parameters are involved in 



Eqs. p2)) and ((T3| . we tested these expressions in various 
instances and they turned out to be always remarkably 
accurate already for i > 5 [13- In Fig. [5] we show C^ and 
R^ for Fq = 0.3, together with a zoom into the long-time 
decay and its comparison with the leading exponential 
decay in the upper inset. The lower inset confirms the 
high quality of the fit of the correction terms to the forms 
given in Eqs. (fT2|) and (fTB]) which are actually indistin- 
guishable from the data. The non-equilibrium coherence 
time T = TTy/T — l/[4arctan(VT — 1)] decreases upon 
increasing |1 — Fo|, i.e., the energy injected into the sys- 
tem and T ^ |1 — FqI^^ for Fq -^ 1. While the parameters 
^i?,C depend on Fg, their ratio Ac/Aji = 1.210(5) does 
not within our numerical accuracy. More details on the 
fitting parameters of the oscillating (lattice) correction 
will be presented in Ref. [19|. 




FIG. 2: (Color online.) Decay of the order-parameter corre- 
lation C^ (blue line) and the linear response R^ (red line) 
for To = 0.3. Upper inset: zoom into the long-f decay that 
demonstrates the exponential relaxation with the characteris- 
tic time r defined in the text (dashed green line). Lower inset; 
(e*''^C+Mc - l)/ac and (e*/^i?^/^fl - l)/ar vs. i; the green 
dashed line is the i~^ " envelope of the damped oscillations, 
in agreement with Eqs. p2p and (|13p . 

The effective temperatures determined in the fre- 
quency and in the time domain arc equivalent in this 
case. As discussed above, replacing /? by a constant /3jfg, 
turns the rhs of the FDT ([S]) in the time domain into a 
series of time derivatives of C^ (i) which yields 



hAR/{2Ac) = tan{hl3^s/2T) 



(14) 



for t — >■ oo, i.e., neglecting the oscillatory corrections 
in Eqs. ([12]) and (US]). Ahcrnatively, Eq. ^ yields 
f]^^(uj = 0) = J^dttR{t)/ J^dtC{t) for a; ^ which 
numerically coincides with the constant value in Eq. (|14p . 
For h(3^g/2T <^ 1 one recovers the classical limit /3^fj ~ 
-R'{t)/[dC'=it)/dt] ~ tAr/Ac. Ah three determina- 
tions of T^ff are shown in Fig. [3] as functions of Fg and 
they are compared to T^ (dashed line) from Eq. ([3|) [7| . 
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FIG. 3; (Color online.) To dependence of the order parame- 
ter effective temperature T^ff compared to TJ^ defined from 
the energy [see Eq. Q] (dashed line). The solid lines, from 
bottom to top, indicate the values determined on the basis 
of the classical limit of the FDR in the time domain, of the 
limit cj — >■ of the frequency-domain FDR, of Eq. H14I) . and of 
the hmit oj — >■ of the frequency-domain FDR but for spins 
separated by a distance r = 10. The dependence on Fq can 
be read from the one of t and Eq. p4|) . 

ogous results, as shown in Fig. [3] for the a;-componcnt 
of two spins at distance r = 10. (Note that differently 
from case r = shown in Fig. ^ correlations with r ^ 
display a light-cone effect due to the finite speed of the 
quasi-particles [i,[l3-[lll-) 

Conclusions. Independently of the functional form of 
the correlations involved, the FDRs allow us to define 
various effective temperatures. We calculated the (self) 
FDR for three observables that are local (cr^'^) or non- 
local (M) in space and local (erf, M) or non-local (trf ) 
in the quasi-particles. af is not compatible with Gibbs 
thermal equilibrium at any effective temperature. The 
frequency-domain FDR for M yields a finite T^g{TQ) in 
the limit oj — > 0. Frequency and time-domain determi- 
nations of T^g{ro) are equivalent. T^ and T^g have the 
same qualitative dependence on Fq but they differ (also 
from T^). This excludes a single temperature effective 
Gibbs description (as the one discussed in Ref. |2J] ) of the 
full stationary dynamics of this model but the question 
remains as to whether some of the temperatures which 
emerge can be attributed a thermodynamic meaning. 

We finally stress that a bona fide thermal behavior 
should be accompanied by the validity of suitable FDTs 
also in the context of quantum quenches. 
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We completed our analysis by studying space- 
dependent correlations and we found that they yield anal- 



[1] A. Polkovnikov, K. Sengupta, A. Silva, M. Vengalattore, 
arXiv: 1007. 5331. J. Dziarmaga, Adv. in Phys. 59, 1063 



(2010). 
[2] D. Karevski, arXiv:cond-mat/0611327. 
[3] M. Rigol, V. Dunjko, V. Yurovsky, and M. Olshanii, 

Phys. Rev. Lett. 98, 050405 (2007). M. Rigol, A. Mura- 

matsu, and M. Olshanii, Phys. Rev. A 74, 053616 (2006). 
[4] M. A. Cazahlla, Phys. Rev. Lett. 97, 156403 (2006); 

A. lucci and M. A. Cazahlla, Phys. Rev. A 80, 063619 

(2009). 
[5] D. Fioretto and G. Mussardo, New J. Phys. 12, 055015 

(2010). 
[6] P. Calabrese, F. Essler, and M. Fagotti, Phys. Rev. Lett. 

106, 227203 (2011). 
[7] D. Rossini, A. Silva, G. Mussardo, and G. E. Santoro, 

Phys. Rev. Lett. 102, 127204 (2009); D. Rossini, S. 

Suzuki, G. Mussardo, G. E. Santoro, and A. Silva, Phys. 

Rev. B 82, 144302 (2010). 
[8] R. Kubo, M. Toda, and N. Hashitume, Nonequilibrium 

statistical mechanics, 2nd ed. (Springer Verlag, 1991). 
[9] L. F. Cugliandolo, J. Kurchan, and L. Peliti, Phys. Rev. 

E 55, 3898 (1997). 
[10] A. Crisanti and F. Ritort, J. Phys. A 36, R181 (2003). L. 

Leuzzi, J. Non-Cryst. Sol. 355, 686 (2009). L. F. Cuglian- 
dolo, arXiv: 1104.4901;L. F. Cughandolo and G. Lozano, 

Phys. Rev. Lett. 80, 4979 (1998); Phys. Rev. B 59, 915 

(1999). 
[11] A. Caso, L. Arrachea, and G. S. Lozano, Phys. Rev. B 

83, 165419 (2011) and refs. therein. 
[12] P. Calabrese and A. Gambassi, J. Phys. A 38, R133 



(2005); F. Corberi, E. Lippiello, and M. Zannetti, J. Stat. 
Mech. P07002 (2007). 

[13] C. Godreche and J.-M. Luck, J. Phys. A 33, 1151 (2000). 

[14] P. Calabrese and A. Gambassi, J. Stat. Mech. P07013 
(2004); J. Stat. Mech. Theor. Exp. POlOOl (2007). 

[15] B. M. McCoy, E. Barouch, and D. B. Abraham, Phys. 
Rev. A 4, 2331 (1971). 

[16] F. Igloi and H. Rieger, Phys. Rev. Lett. 85, 3233 (2000). 

[17] F. Igloi and H. Rieger, Phys. Rev. Lett. 106, 035701 
(2011); arXvi : 1106 . 5248. 

[18] P. Calabrese and J. Cardy, Phys. Rev. Lett. 96, 136801 
(2006); J. Stat. Mech. P06008 (2007). 

[19] L. Foini, L. F. Cugliandolo, and A. Gambassi (unpub- 
lished) . 

[20] C^' (t) has an additional constant term, in and out ol 
equilibrium, equal to (vT + l)~^ for critical quenches. It 
is however irrelevant for the study of FDT. Note also that 
the definition of C± (t) as in Eq. (J4]| requires an extra 
normalization factor L, necessary for global observables. 

[21] T. Niemeijer, Physica 36, 377 (1967). 

[22] P. Deift and X. Zhou, in Singular limits of dispersive 
waves (Lyon, 1991), p. 183, NATO Adv. Sci. Inst. Ser. B 
Phys. 320 (Plenum, New York, 1994). 

[23] The value of r^q for small T quoted sometimes in the 
literature differs from the correct one Teq = 4/(7rT) [23 |. 

[24] A. Mitra and T. Giamarchi, arXiv: 1105.0124. J. Lan- 
caster, T. Giamarchi, and A. Mitra, arXiv: 1105.6039. 



